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ABSTRACT 

We study the effect of rotation on the excitation of internal oscillation modes of a 
star by the external gravitational potential of its companion. Unlike the nonrotating 
case, there are difficulties with the usual mode decomposition for rotating stars 
because of the asymmetry between modes propagating in the direction of rotation and 
those propagating opposite to it. For an eccentric binary system, we derive general 
expressions for the energy transfer AEg and the corresponding angular momentum 
transfer AJg in a periastron passage when there is no initial oscillation present in the 
star. Except when a nearly precise orbital resonance occurs (i.e., the mode frequency 
equals multiple of the orbital frequency) , AEg is very close to the steady-state mode 
energy in the tide in the presence of dissipation. It is shown that stellar rotation can 
change the strength of dynamical tide significantly. In particular, retrograde rotation 
with respect to the orbit increases the energy transfer by bringing lower-order g-modes 
(or f-mode for convective stars), which couple more strongly to the tidal potential, into 
closer resonances with the orbital motion of the companion. 

We apply our general formalism to the problems of tidal capture binary formation 
and the orbital evolution of PSR J0045-7319/B-star binary. Stellar rotation changes 
the critical impact parameter for binary capture. Although the enhancement (by 
retrograde rotation) in the capture cross section is at most ~ 20%, the probability 
that the captured system survives disruption/merging and therefore becomes a binary 
can be significantly larger. It is found that in order to explain the observed rapid 
orbital decay of the PSR J0045-7319 binary system, retrograde rotation in the B-star 
is required. 



Subject headings: binaries: close - pulsars: individual (PSR J0045-7319) - stars: 
neutron - stars: oscillations - stars: rotation - hydrodynamics 
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1. Introduction 

The standard weak friction theory for tidal interaction in binary stars, first introduced by 
Darwin (1879), and developed in detail by many authors (e.g., Alexander 1973; Kopal 1978; Hut 
1981), is based on the assumption of static tide in hydrostatic equilibrium. In the presence of 
finite dissipation, a tidal lag develops between the tidal biilge and the external tidal potential, 
and the resulting tidal torque drives the spin and orbital evolution. This theory has been applied 
successfully to planet-satellite systems (e.g., Goldreich &: Peale 1968), and can also be used to 
describe binaries of low-mass (late-type) stars if appropriate turbulent viscosity due to convective 
eddies is used (Zahn 1977; Goodman & Oh 1997). The static tide approximation, however, is not 
appropriate for many other binary processes, especially those involving massive (early-type) stars 
with radiative envelopes, where dynamical excitation and radiative damping of low-frequency 
g-modes dominate the tidal interaction (Zahn 1977; Goldreich & Nicholson 1989). In these early 
works, the effect of rotation on the dynamical tide is neglected or introduced phenomenogically 
(but see Nicholson 1979 and Rocca 1989, where some aspects of rotational effects are explored), 
and only circular (or near circular) binary orbits are considered. 

In this paper we study dynamical tides in rotating, eccentric binary stars. Our work is 
motivated by the recent observations of the PSR J0045-7319/B-star binary (Kaspi et al. 1996), 
which contains a radio pulsar and a massive rapidly-rotating B-star in an eccentric 51 days orbit, 
with periastron separation of only 4 stellar radii. This system therefore constitutes an excellent 
laboratory for studying dynamical tidal interaction. In a preliminary analysis (Lai 1996) of the 
rapid orbital decay of the binary as revealed by timing observation (Kaspi et al. 1996), it was 
realized that rotation can have rather significant effect on the strength of dynamical tides. In 
particular, when the B-star has a retrograde rotation (with magnitude of one half of the maximum 
value, for example) with respect to the orbital motion, the dynamical tidal energy is increased 
over the nonrotating value by two orders of magnitude. Roughly speaking, this comes about 
because retrograde rotation shifts the tidal resonance (where the mode frequency equals the 
external driving frequency) from high-order g-modes to lower-order ones, which couple much more 
strongly to the tidal potential. It is this enhancement of dynamical tide by retrograde rotation 
that allows the binary orbit to decay on a short time scale of 0.5 Myr, although it has been 
suggested (correctly in our opinion) by Kumar and Quataert (1997; hereafter KQ) that differential 
rotation is also required in order to keep the mode damping time relatively short (see §6 for 
more discussion). In this paper we present a full analysis of dynamical tides in rotating stars, 
with applications to the orbital evolution of the PSR J0045-7319/B-star binary. We shall focus 
on the d3mamical aspects of the problem, i.e., those aspects that are independent of dissipation 
mechanisms. We particularly emphasize the difficulty with the usual mode decomposition because 
of the asymmetry between modes with opposite signs of pattern speed. Such difficulty has led to 
quantitatively misleading results concerning the tidal energy in stars with rapid prograde rotation 
(KQ; see §6.3 for discussion). 



Another prominent example in which dynamical tides play an important role is binary 
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formation from tidal capture. This process was originally suggested by Fabian, Pringle &; Rees 
(1975) to explain the origin of low-mass X-ray binaries in globular clusters: During a close passage 
of two stars in an unbound orbit, the excess orbital energy is transferred through dynamical 
tidal excitation to the internal stellar oscillations, resulting in a bound system. Although it has 
been recognized in recent years that primordial binaries play an important role in the dynamics 
of globular clusters, and may also be involved (via exchange interaction with another star) to 
produce low-mass X-ray binaries and other stellar exotica (e.g.. Hut et al. 1992; Davies 1996), 
tidal capture can still be an efficient means to produce tight binaries, thereby providing an energy 
source to reverse core collapse. Many theoretical issues related to tidal capture process have been 
studied, including linear and nonlinear energy transfer (Press and Teukolsky 1977, hereafter FT; 
Lee & Ostriker 1986; Khokhlov, Novikov k, Fethick 1993), evolution of the binary subsequent to 
capture (Kochanek 1992; Novikov, Fethick & Folnarev 1992; Mardling 1995) and dissipation of 
tidally excited modes (Kumar &; Goodman 1996). However, the possible role of stellar rotation on 
the tidal energy and angular momentum transfer has not been considered. Indeed, we find in this 
paper that rotation changes the cross-section for tidal capture, and can significantly increase the 
probability of forming binaries through tidal interaction. 

Our paper is organized as follows. In §2 we present the basic equations describing the tidal 
excitations of normal modes in a rotating star and the orbital response to the excited oscillations. 
We then derive in §3 general expressions for energy transfer and angular momentum transfer in 
a periastron passage when there is no initial oscillation. For an eccentric system, this energy 
transfer is directly proportional to the steady-state mode energy in the presence of dissipation, 
and provide a measure of the strength of the dynamical tide (see §6.1). In §4 we examine the 
properties of rotating modes that are relevant for tidal excitation; these are used in specific 
calculations presented in §§5-6. We apply our general formalism to the tidal capture problem in 
§5 and the orbital evolution of FSR J0045-7319/B-star binary in §6. A brief discussion of related 
problems is given in §7. 

2. Basic Equations 

Consider a rotating star with mass M, radius R and spin fig, interacting with a companion 

M' . We treat M' as a point mass, although generalization to the case where both stars have finite 
structure is straightforward. For simplicity here we assume the spin axis is perpendicular to the 
orbital plane (but Vlg can have both positive and negative signs, corresponding to prograde and 
retrograde rotations with respect to the orbit). Generalized equations for arbitrary spin-orbit 
inclination angles are summarized in Appendix A. 

In the linear regime, the perturbation of the tidal potential on M is specified by the 
Lagrangian displacement £,{v,t) of a fluid element from its unperturbed position. In the inertial 
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frame, the equation of motion can be written in the form: |^ 

gf+2(v.V)§ + C.f=-VC. (1) 

where v = $7^ x r is the unperturbed fluid velocity, C is a self-adjoint operator (see Lyden-bell &; 
Ostriker 1967). The external potential in Eq. (||) is given by 

CM' 

U{r,t) = - p = -GM'J2Wim^exp{-im'^)YU9,4>), (2) 

where we have chosen a spherical coordinate system centered at M with the z-axis long the 
rotation axis and the x-axis pointing to the periastron: r is the position vector of a fluid element 
the star, and D(t) = [D{t),7r/2,^{t)] specifies the position of the point mass M' [D{t) is the 
orbital separation, ^{t) is the orbital phase], and Wim is a coefficient as defined in PT: 

1/2 



(3) 



[here the symbol (— is to be interpreted as zero when k is not an integer]. In equation (|2|), 
the I = and I = 1 terms should be dropped, since they are not relevant to tidal deformation. 
For 1 = 2, the only nonzero W^m's are W2±2 = (37r/10)^/2 and W20 = reflecting the 

symmetry of the tidal potential with respect to (p ^ (p + n. Thus only modes with m = 0, it 2 are 
excited. 

An eigenmode ^Q,(r, i) = ^Q(r)e*'^"* oc e*°""*+*'^'^ satisfies the equation 

- alic. + 2i(T„ (v • V)^"; + C • = 0, (4) 

where aa is the mode angular frequency in the inertial frame, and a = {njm} specifies the mode 
index: n gives the number of nodes in the radial eigenfunction (the order of the mode) , j specifies 
the number of nodes in the 0-eigenfunction (j reduces to I in the nonrotating limit), and m is the 
azimuthal mode index. We shall use the convention > in this paper, so that a m > mode 
has a retrograde pattern speed {—a/m < 0), while a m < mode has a prograde pattern speed. 
Although the normal modes are not complete in strict mathematical sense, any reasonable initial 
data for Eq. (|^) evolve in time as a linear superposition of the eigenmodes ^a(r) (see Dyson & 
Schutz 1979 for an analysis of the completeness problem). Thus we write C(r,i) = J2a(^a{t)^a{Y), 
and normalize the eigenmode via Jd^x p^a • Co = 1- Substituting this expansion into Eq. (||) and 
using Eq. (^, we obtain the dynamical equation for the mode amplitude: 

■■i2 o-D/- • \io- GM'WimQal -im* /r;\ 
Ua + a^aa - 2tBa{aa - tdaaa) + = 2^ J^l^ 6 , (5) 



^The mode amplitude can be equivalently studied in the rotating frame. In this case, the equation of motion is 
d'^i/dt^ + 2Q.S X (d^/dt) + C' • ^ = -V!7, and in Eq. (|) we should replace $ by (<1> - Q.st). Our treatment in the 
inertial frame is more convenient for generalizatiion to the case of arbitrary spin-orbit inclination angles (Appendix 
A). 
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where Qai is the tidal couphng coefficient (the overlap integral) defined by 

Q^i = J d^xpC ■ V(rVz^) = J Sx5pl{rWira), (6) 

[here 5pa = —V • (p^a) is the Eulerian density perturbation], and the function is given by 

Ba, = ij d^xpC ■ (v • V)^"; = -rnQs + i J d^xp^ ■ i^s x (7) 

Note that p^^ • (v • V).^^ = V • (/Ov||^p) — |^opV • (/ov) — p^a • (v • V).^*, and one can easily show 
that Ba is real. In deriving Eq. (||), we have assumed that different modes are orthogonal to each 
other, and have neglected the off-diagonal terms Ba/s = i J d?x p^^ • (v • V)^q. This approximation 
allows us to decompose ^(r, i) into normal modes and obtain a simple equation for individual 
mode amplitude. The limitation of this approximation will be examined in the next section [see 
the discussion following Eq. (28)]. 

We have also introduced in Eq. (^) the quantity 7^, which is proportional to the mode 
damping rate. Multiply Eq. ^ by a* and take the real part, we find: 

-^--27„|aa| +2^ a^e , (8) 

where the mode energy is given by |^ 

= ^|aa|2 + ^((j2 -2S„a„)|a„|2. (9) 

The last term in Eq. (^) gives the rate at which energy is transferred to each mode. For free 
oscillation (without external driving), oc e*'^"*, the mode energy becomes cra{(^a — Ba)\aa\'^, and 
Eq. (^) reduces to dEa/dt = —2TaEa, with the mode amplitude damping rate given by 

r„ = (10) 

The excited modes affect the orbital motion through the tidal interaction potential 

Vu.e = -GM' jd^x'^ = - E^^^^5^«^We^-*. (11) 
Thus the equations of motion for the orbit are 

al 

d_{p^ = Y^^m^M:^^e^--a.[t). (13) 

al 



^The canonical energy associated with perturbation ^(r, t) is generally given by Ec = (1/2) Jd^s p ([S^/c^fp +5* 
C ■ ^). For i) = aa{t)£,a{v), the canonical energy reduces to Eq. (H). 
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Given the mode properties {aa,Qai), equations (|5|),([l^ and ( |T^ ) completely determine the 
time evolution of the tidally excited oscillations and the binary orbit. In general, the energy 
transfer between the star and the orbit in a periastron passage depends on the phases of the 
oscillation modes and thus varies from one orbit to another (e.g., Kochanek 1992; Mardling 1995). 
For a highly eccentric system (such as parabolic tidal capture and the PSR J0045-7319 binary, to 
be discussed in §5 and §6, respectively), an important measure of the strength of dynamical tide 
is the energy transfer to the star, AEg, and the correspinding angular momentum transfer, AJg, 
during the "first" peraistron when there is no initial oscillation. General expressions for AEs and 
AJs are derived in the next section (§3). The relation between AEg and the steady-state tidal 
energy (Eg) in elliptical binary system in the presence of dissipation is established in §6.1. 



3. Energy Transfer and Angular Momentum transfer 

Here we derive the expressions for the energy transfer AEg and angular momentum transfer 
AJg to the oscillation modes in a single periastron passage assuming the star has no initial 
oscillation. The binary orbit is characterized by the periastron distance Dp and two dimensionless 
quantities: the eccentricity e, and the ratio between the periastron passage time and the dynamical 
time of the star: 

where Mt = M + M' is the total binary mass, and Op is the orbital angular frequency at periastron 
in units of {GM / B?)^/"^ . We outline two derivations of AEg and AJg. 

Our first derivation is patterned after PT. Since for very eccentric system, the orbital 
period is much longer than the timescale for periastron passage, and thus much longer than the 
periods of the tidally-excited normal modes, we can consider a continuous Fourier transform, viz. 
aa{t) = Jda e^'^^aa{(^)- Equation (|5|) then becomes 

GM' 

i^^l - + 2Baicr - (7a)da + 2i7„<Ta„ = —j^QalKimicr), (15) 

I 

where Kimic) depends on the orbital trajectory: 

Ki^(a) = / dt(-^] exp(im$ + iat), (16) 



2tt J-oo \D 

with the integration centered around the periastron (For elliptical orbit, the upper limit and lower 
limit of the integration should be replaced by Porb/2 and — Porb/2). The mode amplitude is then 
given by 

a.{t) = y Idae^'^'^, GM'Q^^KU^) 
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The energy transfer AEg to the oscillation modes (in the inertial frame) is obtained via 

AE, = -Jdt Jd'xpf^-VU*= Jdtj:^^^^^^e^-^d^{t), (18) 

al 

(only the real part of the expression is meaningful). Using Eq. (|T7|), this becomes 

, ^ V- {GM'fQc,lQal'WlmWl>„^ f f , eim${t)+im<I.(t') 



all 

where 

A^{t,t')= da- , — . (20) 

This integral can be easily evaluated using a countour in the complex plane, yielding 

A^{t, t') = ^ [a„e^""(*+*') + {aa - 2S^)e~^('^~2B.){t+t')l ^ (21) 

for t + t' > 0, and Aa{t, t') = for t + t' < 0. The energy transfer is therefore 

Air. = — i:(-j T„,(,.e), (22) 

with 



7]r(f?, e) = vr^ V — ^ [(TQ-R'im(o-Q)-f^«'m(o-a) + (cTa " 2Ba)Ki_rniaa " 2Ba)Ki, _rn{(^ a " 2-B„)] . 

(23) 

In writing Eq. (p^), we have factored out all dimensional physical quantities so that Tui is 
dimensionless, i.e., all quantities in Tui are in units such that G = M = R = 1. The angular 
momentum transfer is obtained via 

AJg = Jdt J d'x6p (-^) = / dt^{-^m f^'^lf'^U ^-^a^{t). (24) 
Using a similar procedure we find: 



with 



Swiv, e) = vr^ V ^„ [Er/m(c^a)i^i'm(c^a) " ^^^-^.(o-a " 2Ba)Kii^m{aa - 2B„)] . (26) 

a <^a — Ba 

For a given a = {njm}, the first terms inside the square brackets of Eq. ( p3| ) and Eq. ( p6[ ) 
correspond to a mode with frequency iTq, while the second terms correspond to a mode with 
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frequency — ((Jq — 2Ba). To see this more clearly, consider another derivation of AEg and AJg. 
We simply integrate Eq. ^ directly to obtain aa{t): 



^ -im^(t')-i(7at' rt ^-im^(t')+i(cra~2Ba)t' 



. . ^ ^ GM'WlmQal 

(27) 

where we have neglected the damping term by setting Fq, = 0. Suppose the periastron passage 
occurs at t = 0, then for t >> 0, we have 

= E ■J^u^'^'^' ^ U^'-'KU'Tc) - e-^(--2^'^)*i^z_^(cT, - 2i?„)l . (28) 
^ iDl+\aa- B^) L J 

Using Eq. (|9|) we find that the total mode energy after the periastron passage is again given by 
Eqs. (p^ 



Equation ( |28| ) clearly reveals that after excitation in a periastron passage, a mode ^{r, 9)e'^"^'^ 
oscillates at frequencies CTq, (its eigenfrequency) and —{a a — 2Ba). How can a mode oscillate at 
a frequency other than its own eigenfrequency? This conundrum comes from our approximate 
treatment of mode decomposition leading to Eq. (|5|). Note that the mode ^q, oc e~*°"*~*™''^ is 
physically identical to the mode q^^+i-^'P ^ Yoi a given n,j and \m\ = 2 (for example), there 
are four mathematically distinct modes: | 

imir^t) = ^2{r,e)e'^'^+"'^\ {r,e)e-''^'^-"'^\ (retrograde mode); (29) 

i^{Y,t) = f_2(r,^)e-*2'^+*'"-2*, |l'2(r,0)e^2<^-^'"-2*, (prograde mode), (30) 

and only two of them are physically distinct, one with prograde pattern speed and another with 
retrograde pattern speed. Clearly, within the context of our approximation, we should identify 
(dm — 2Bm) with U-m (where the other mode indices nj are surpressed). This (— (T_m) mode has 
wavefunction ^^^{r,6), different from ^m{r,9), and thus its ampitude should be proportional to 
Q-m rather than Qm as indicated in Eq. (|2^). But this is an artifact of our approximation. In the 
slow rotation limit {\i^s\ « CTa), the mode frequencies are related to the nonrotating value a^'^^ 
by dm = cr^'^^ + Bm and cj_m = cr^'^^ + B-m = cr^^^ — Bm, while the wavefunctions are unchanged 
to the leading order (so that Qm = Q^m = Q^°^)- Thus Eqs. (|2|), (||), (|2|), (|2|) and (|2|) are 
exact in this limit. For larger rotation rates, they should still be good approximations as long as 
am — 2Bm is close to cT_m, and the difference between Qm and Q-m is small — We find that these 
are satisfied for the relevant modes considered in §5 and §6. 

The preceding discussion leads us to adopt the following ansatz: The amplitude of each mode 
(specified by a = {njm} and the direction of pattern speed) is given by the first term of Eq. (p7|), 
while the contribution of the second term to AJs and AEg is identical to that of a corresponding 



^It can be shown that if f(r, 6')e""'^+"'* is a solution of Eq. (|), then £^'{r, 
^' — (^,.,Cfl) ~^4i)- These two modes have the same 5p(r,9). See Schutz (1979). 
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(— m) mode with the same pattern speed. The dimensionless energy trnasfer T and angular 
momentum transfer S are therefore given by 

Tii'{'n,e) = 27r^V (7aKimiaa)Ki>miaa), (31) 

Sw{v,e) = 27r^ ^ ^"'^^^^'"^^ KUaM'^CTa), , (32) 

where the factor of 2 results from the equivalence of the (m, cJq,) mode and the (— m, — cTq,) mode. 
Since the mode pattern speed is {—aa/m), the contribution of each mode to the energy transfer 
AEa and angular momentum transfer AJq are related by AJq = {—m/ aa)lS.Ea, as seen from 
Eqs. (|3T|) and (p^). It is also useful to consider the total mode energy /S.E^'^ in the star's rotating 



frame. Since the rotational kinetic energy associated with the angular momentum transfer is 
rj^AJs, we have IS.E^^ = AEg — fi^AJs, and the corresponding dimensionless energy is 

^iPiv, e) = Tii'{r], e) - ^sSwirj, e) = li? V '^"'^p' uJaKimicra)Ki>„,{aa), (33) 

a - Ba 

where cj^ = (Tq + mi^s is the mode frequency in the rotating frame. We emphasize that, although 
our procedure leading to Eqs. (|3l|)-(33) is somewhat ad-hoc mathematically [The "rigorous" 



derivation of Eqs. (23), ( p^ and ( [2^ ) are exact only in the slow-rotation limit], the forms of these 
expressions are naturally expected on the physical ground: the retrograde mode should be related 
to |Q2-f^i2(c2)P and the prograde mode to |(5-2-^/-2(f'"-2)P- The only uncertainty is the factor 
(cTq, — Ba)/iOa = 1 — (i/uJa) J d^x • (fig X ^a) [cf. Eq. (^)] which only serves as a insignificant 
"correction" to Qa (see §6 for further discussion). 

Although the above expressions apply for general /, in this paper we shall consider only the 
quadrupolar {I = 2) tides. The contribution of higher-/ tidal potential to AEg is of order [R/Dp)^ 
smaller than that of the leading I = 2 term. 



4. Properties of Modes in Rotating Stars 

We consider two simple stellar models in this paper: (i) A F = Fi = 5/3 polytrope (T is the 
polytropic index, Fi is the adiabatic index), representing low-mass (M ^ 0.5Mq) MS star with 
thick convective envelope; (ii) A F = 4/3, Fi = 5/3 polytrope, representing high-mass (larger than 
a few solar mass) star for which the envelope is radiative; it is also a good approximation to the 
outer structure of intermediate-mass stars (~ \Mq). 

We examine the effects of rotation on the mode properties (frequency and tidal coupling 
coefficient) relevant to tidal excitation. 
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4.1. f-Mode and p-Modes 

For low- mass, convective star (which does not have any g-mode), the f-mode and lowest order 

p-modes are the most important modes that absorb tidal energy. For Clg much smaller than the 

mode frequency, perturbation theory is adequate (e.g., Unno et al. 1989). To leading order, the 

mode wavefunction is unchanged, and the mode frequency is modified from the nonrotating value 
(0) u 

<ynlm = (^ni + Bnlm = (^nl ~ "^^s + mCnl^s, (34) 

where rnQgCnl = ^ Jd^^P^a ' ^ Using the zeroth order cigcnfunction 

Lmir) = [Ur)er + ^±{r)V±\Yim{e,4>), where V = egid/dO) + (e<^/sin^)(5/9(/)), we 
have ^ 

Cni= I pr\2iri^ + e±)dr, (35) 
Jo 

and the eigenfunction is normalized via Jd^xp^* ■ ^ = Jdr pr'^[^^. + l{l + 1)^^] = 1. In the 
rotating frame the mode angular frequency is tOnim = ^nl + 'f^^nl^s- For the F = 5/3 polytrope, 
we have C02 = 0.4955 (f mode), C\2 = 0.1525 (pi mode), and C22 = 0.0787 (p2 mode); for 
the F = 4/3, Fi = 5/3 polytrope, we have C02 = 0.2543 (f mode), Cu = 0.1538 (pi mode), 
and C22 = 0.0818 (p2 mode). For higher Qg, linear perturbation theory breaks down, but 
accuarte numerical results are available (e.g., Managan 1986; Ipser &; Lindblom 1990. The later 
is essentially an exact method). In our calculation of tidal capture below (§5), wc shall adopt the 
linear approximation at small 0,s with correction for high 0,s using the results given in Managan 
(1986). The modification to the mode cigcnfunction is negligible and therefore we set Qa = Qa2 
to be the same as that of a nonrotating star. 



4.2. g-Modes 

For intermediate to massive stars with radiative envelopes, g-modes dominate the tidal energy 
transfer because their relatively low frequencies match the orbital frequency at periastron (see 
§§5-6). Since we shall consider the cases where the rotation frequency is comparable or larger than 
the g-mode frequency, the perturbation theory is not adequate. An approximate treatment of 
g-modes in rapidly rotating star is based on the so-called "traditional approximation" (Chapman 
& Lindzen 1970; Unno et al. 1989; Bildsten et al. 1996), where the centrifugal distortion of the 
star is neglected, as well as the Coriolis forces associated with the horizontal component of the 
spin angular velocity. Note that this approximation is strictly valid only for high-order g-modes 
for which Ua « N (where N is the Brunt- Vaisala frequency) and >> ^r, and for « ^ 
(so that the radial component of the Coriolis force can be neglected compared to the buoyancy 
force), but it also provides a good estimate even for low-order modes. 

Neglecting the perturbation in the gravitational potential (Cowling approximation) and 
adopting the traditional approximation, the radial Lagrangian displacement and Eulerian pressure 
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perturbation can be written as 



rn(j> 



(36) 



where Hjm{Q) is the Hough function, satisfying the Laplace tidal equation (Chapman & Lindzen 
1970). Separating out the angular dependence, the fluid continuity equation and Euler equation 
(in the rotating frame) reduce to a set of coupled radial equations: 



d 
dr 



-^5P + p{uj^ - N^)ir, 



cj cjp uj^p 



5P, 



(37) 
(38) 



where the mode index a = {njm} has been surpressed, and (7 > is the gravitational acceleration, 
Cs is the sound speed. The eigenvalue of the Laplace tidal equation, A, depends on m and the 
ratio q = 2^}s/^a- For q (nonrotating case), the function Hjm{0)e^"^'^ becomes Yim{0,4>) while 
A degenerates into /(/ + 1). 

The properties of the Hough function have been extensively studied (Longuet-Higgins 1967). 
We have adopted the numerical approach of Bildsten et al. (1996) in our calculation. For given 
q and m, the eigenvalue A is obtained by solving the angular equation. The radial equations 
(13)" (H) then solved together with appropriate boundary conditions to obtain uj, from which 
the actual rotation rate fls = q'-o/2 is recovered. Figure 1 shows the frequencies (in the rotating 
frame) of several j = 2 g-modes against the rotation rate (Both uj and 0,s are plotted in units 
of uj^^\ the mode frequency at zero rotation). These modes have ? = 2 in the 0^ = limit. 
For high order g-modes, a WKB analysis for the radial equations gives uj oc -v/A. Thus we have 
uj/u) 



A/6 and ^Is/^ = qV A/24. For | Os/ct; | < 1.1 (which is satisfied by g-modes with 



(0) 



n < 13 when < 0.5), our numerical results for the g-modes of the F 
can be fitted by the following analytic expressions to within 1%: ^ 



4/3, Fi = 5/3 polytrope 



UJ 



1 + 

0) = 1 + 

= 1 + 



(l/3)f]s + (13/42)0^ - 0.064J^^-^ /, 



{6/7)ni-o.3i\nsr'\ f, 

-(1/3)17, + (13/42)f^2 -O.llSf^^l /, 



(m = 2), 
(m = 0), 



[m 



-2), 



(39) 
(40) 
(41) 



where a) = uj/uj^'^\ D,s = i}s/^^^\ and the factor / < 1 depends on specific g-modes: In the 
WKB limit (high-order g-modes with n ^ 00) we have / = 1. The values of / for other 
selected modes are: / = 0.48, 0.81, 0.9, 0.95 for gi, g^, gio, g2o- Equations (^)-(|4l|) break 
down when ^ 1.1. In this high-Os regime, the following asymptotic expressions can be 
applied: O ~ 1.54-\/n7 (for m = 0, 2), and uj ~ 0.82 (for m = —2). Note that these expressions 



^ In these expressions, the Unear and quadratic terms are derived from an analytic expansion of A for small q 
(Bildsten & Ushomirsky 1996, private communication), and the last terms are based on numerical fitting. 
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[Eqs. (^)-(|4l|)] are valid only for Qg > 0. The eigenfrequencies for < can be obtained using 
the relation ujm{^s) = ^-mi—^s)- 

Rotation also changes the tidal coupling coefficient Qa2, although the correction is not 
significant since all the g-modes we include in our calculation (n < 20) that are strongly excited 
satisfy 0,s/^^^^ < 1-7, i.e., 0^ is not much larger than the mode frequency. Expressions for 
evaluating the Qai and Bq, are given in Appendix B. Figure 2 shows the coupling coefficient 
Qa = Qa2 for several g-modes. 

For rotating stars, the j = 4, 6, 8, • • • (corresponding to / = 4, 6, 8, • • • in the ^ limit) 
modes with m = 0, ±2 are also coupled to the quardrupolar tidal potential. (The coupling with 
the j = 3, 5, 7, • • • modes vanishes because the eigenfunctions are odd with respect to the equator 
while the tidal potential is even). The coupling goes to zero at 0,s = 0. Some properties of the 
j = 4 modes are plotted in Figure 3 and compared with the j = 2 modes. Two factors make these 
higher-j modes less important to energy transfer during a tidal encounter than the j = 2 modes: 
(i) The J > 4 modes have higher frequencies than the j = 2 modes of the same radial order, 
therefore resonance condition (see next section) is satisfied for higher-order modes, which couple 
less strongly to the tidal potential; (ii) The j > 4 modes have much smaller coupling coefficients 
than the j = 2 modes of the same radial order (The angular integral Qg is plotted in Fig. 3(b); 
the radial integral Qr is slightly smaller for the j = 4 modes than for the j = 2 modes). We have 
checked that the contribution of the j > 2 modes to the energy transfer is always less than a few 
percent. Thus they are neglected in our calculations below. 

4.3. r-Modes 

For nonrotating stars, these are "trivial" toriodal modes with zero frequency and zero density 
perturbation, and thus they do not couple to the tidal potential. With rotation these modes are 
mixed with spheroidal components, and attain finite frequencies. The r-mode (sometimes also 
called "quasi-toroidal mode") is the generalized form of Rossby waves on a spherical surface as 
studied in terrestrial meteorology. For small rotation rate \0,s\ « !> r- modes can be studied using 
perturbation expansion in powers of ilg (Provost et al. 1981; Saio 1982). To linear order in fi^, the 
mode frequency increases as a; = 2mQs/l{l + 1), while its Lagrangian displacement remains purely 
toroidal. Finite density perturbation comes in only in the second order of Vlg. For higher O^, 
r-modes can also be associated with eigenvalues of the Laplace tidal equation (Longuet-Higgins 
1967). Since the coupling between r-modes and the tidal potential is weak {Qa oc fi^), we shall 
neglect them in our calculations below. 
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5. Tidal Capture Binary Formation 

We now apply the results of §§3-4 to study the effect of rotation on tidal capture. As long as 
the relative velocity between the stars at infinite separation is much less than the escape velocity 
from the star, a parabolic orbit is a good approximation. We consider two representative stellar 
models in §5.1 and §5.2. 

5.1. Intermediate to Massive Stars: F = 4/3 Polytrope 

For the F = 4/3, Fi = 5/3 polytrope, the functions T^''^ = T§ , S2 = S22 and T2 = T^''^ + ^1^82 
[where 0,s is the rotation rate in units of (GM/ii^)^/^] are plotted against rj in Figure 4 for 
^Is = 0, ±0.2 and ±0.4. We clearly see that for < (retrograde rotation), the energy and 
angular momentum transfers can be significantly increased over the nonrotating values. This 
comes about for the following reasons: During a perastron passage, the most strongly excited 
modes are those (i) propagating in the same direction as the orbital motion of the companion 
(corresponding to the m = —2 modes in our notation), (ii) having frequencies in the inertial frame 
comparable to the "driving frequency", which is equal to twice of orbital frequency at periastron, 
i.e., (Jq = oJq — mQs = 2^2^ (the "resonant condition"), and (iii) having relatively large Qa- Since 
higher-order (lower frequency) g-modes have smaller couping coefficients than the low-order ones, 
the trade-off between (ii) and (iii) implies that the modes that absorb most tidal energy are 
those with frequencies higher than the resonant mode. For example, when 77 = 7 and fig = 0) 
the resonant mode is gir^, with frequency cua — 0.4 ^ 2Qp, while the dominant modes in energy 
transfer are gg-gs (which have Ua = 0.85 — 0.68). A retrograde rotation "drags" such m = —2 
waves backward, so that the mode frequencies in the inertial frame arc lowered. As a result, the 
dominant modes in energy transfer are shifted to lower radial orders. Because Qa increases rapidly 
with decreasing mode order or increasing mode frequency (we find from numerical calculations 
that Qa oc for nonrotating F = 4/3 polytrope), the energy transfer is greatly increased. For 
example, at = —0.4 and rj = 7, the dominant modes are g3-g5, the net energy transfer (in the 
inertial frame) AEg is about 20 times larger than the nonrotating value. 

While relatively small prograde rotation (Jl^ > 0) decreases the energy transfer, we see from 
Fig. 4 that for sufficiently large 0,s (or more precisely for sufficiently large ratio of Qg/Qp), AJg 
or even AEg can become negative, i.e., angular momentum and energy arc transferred from the 
star to the orbit (see also Fig. 8 in §6). This behavior can be understood as follows: The angular 
momentum of a mode is related to its energy in the rotating frame by Ja = {—m/oJa)Ea \ thus the 
prograde mode (m = — 2) carries positive angular momentum, while the retrograde mode carries 
negative angular momentum (Note that Ea"^ is always positive). For a fixed t], prograde rotation 
increases the frequencies of the prograde modes, shifts the resonance to higher radial orders (which 
have smaller Qa), therefore decreases the positive angular momentum transfer to the star. On the 
other hand, the retrograde modes are "dragged" forward by the rotation, and their cr's become 
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smaller and may even become negative (i.e., they are prograde in the inertial frame). Therefore 
their contributions to the energy transfer become increasingly more important and the negative 
angular momentum transfer to the star increases as increases. When the negative angular 
momentum transferred to the m = 2 modes becomes larger than the positive angular momentum 

(r) 

transferred to the m = —2 modes, the net AJg changes sign, and at the same time, AEs begins 
to increase with increasing (for a given rj). At even larger 0^, when the change in the rotational 
energy ^IgAJs < (due to angular momentum lost to the orbit) dominates over the kinetic energy 
in the modes, the net AEg becomes negative. Our numerical results for the critical fi^ at which 
AJs or AEs changes sign can be fitted nicely by (for 0.1 < J^p < 0.35) 

= 1.2J^p + 0.08, for AJs = (42) 
Qs = 1.51^p + 0.07, for A^, = 0. (43) 

Kumar &: Quataert (1997) claimed that AJs and AEs become negative when il^ ^ rip. This 
is different from our Eqs. (|4^)-(|4^. We believe this quantitative difference arises from their 
problematic treatment of mode decomposition. We defer a discussion of this point to the next 
section (§6.3). 

The maximum Dp for capture is determined by the condition AE^^\M) + AE^^\M') = fiv'^/2, 
where AE^\M) = AE^'^ and AE^^\M') are the energy transfer to M and M' respectively, ^ 
is the reduced mass and v is the relative velocity of the two stars at infinite separation. For 
1.8 < T] < 3.5 (the relevant parameter regime for capture) the function T2 can be fitted to the form 

T2 = 4''^ + nsS2 = Ar]-'', (44) 

where the fitting parameters A, a for several 0,s are listed in Table 1. Let q = M' /M and 
A = 1 + AE^'^\M') / AE^^\M) , we then obtain the critical periastron distance for capture: 



^cap — ^ 



2Aylg(l + 



U2+q)/2 



2/3(4+a) 


V 







-4/3{4+a) 

(45) 



Figure 5 depicts numerical results for the encounter between a IMq star (modeled as a 
r = 4/3, Fi = 5/3 polytrope) and a M' = I.4M0 point mass (neutron star). With A = 1, 
Equation (||) and Table 1 give Dc^^/R = 2A2v^^ '^^^ for fi^ = [This agrees with Lee & Ostriker 
(1986)], L>cap/^ = 2.10 wj-Q^-^^^ for ^s = 0.5, and D^^^/R = 2.77 v^^-'^^^ for = -0.5, where 
v\Q = t;/(10kms^-'^). 



5.2. Low Mass stars: F = 5/3 Polytrope 

This is more applicable for low-mass (convective) stars (e.g., those in globular clusters). In 
this case, only the f-mode (u; = 1.456 for ^Ig = 0) is strongly excited, while the contribution of 
p-modes to the energy transfer is always less than a few percent because of their larger frequencies 
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and smaller Qa (e.g., for = Oj the pi mode has u = 3.21, implying that it is almost always 

out of resonance, and its Qq. is an order of magnitide smaller than that of the f-mode). We find 
similar rotational effect on energy transfer as in §5.1, although it is less prominent (since there is 
only one radial mode involved). For sufficiently large 0^, the angular momentum transfer AJg 
changes signs, but this occurs only at Qg/^p ^ 3.5; the energy transfer AEs never changes sign as 
long as we restrict to the regime fig < 0.6 (i.e., fig is less than the maximum rotation rate). 

Motivated by the functional form of Kij^ in the asymptotic limit [rju >> 1; see Appendix C), 
we can fit the function T2 to the form: 

T2 = T^'^ + 0,52 = Arf exp(-a77), (46) 

with the fitting parameters A and a given in Table 1. This fitting is valid for the range of r/ near 
tidal capture. In the asymptotic limit, including only the m = —2 f-mode contribution, we have 
a = 4v^o"/3. For Q.s = 0, this gives a = 2.74. Figure 6 shows numerical results for encounter 
between a IMq star (modeled as a F = 5/3 polytrope) and a M' = 1.4M0 point mass (neutron 
star). Again, we see that the critical capture radius increases with increasing retrograde rotation. 

5.3. Discussion on Tidal Capture Binary Formation 

The results of §§5.1-5.2 indicate that stellar rotation can increase or decrease the critical 
capture radius by up to ~ 20% (for nearly maximumly rotating stars), leading to similar change 
in the capture cross-section (Jcap — 27rG(M -|- M')Dp/v'^ (where v is the relative velocity of the 
stars in infinite separation). In globular clusters, the stellar rotation is likely to be much less than 
the maximum since the stars persumably have slowed down significantly by magnetic breaking. 
For open clusters (with much lower velocity dispersion, ~ 1 km s~^) or other young star clusters, 
such rapid rotation is certainly realistic (e.g., the typical surface rotation velocity for B stars is 
~ 400 km/s, corresponding to a rotation rate close to the maximum). While the 20% change in 
-Dcap may seem insignificant by itself, we note that the enhancement to the probability for tidal 
capture binary formation can be much larger. Indeed, it has been suggested that almost all tidally 
captured bound systems may merge during the subsequent periastron passages, leaving no binaries 
behind (e.g., McMillan et al. 1987; Kochanek 1992). The critical periastron distance for merging, 
-Dmergei is Uncertain, but a reasonable lower limit is given by the sum of the radii of the two 
stars. Clearly, if Dmerge is close to -Dcapj then a relatively small increase in £)cap can lead to much 
larger increase in the efficiency for binary formation. If so, we can predict that binaries formed by 
tidal capture must perferentially have retrograde rotations. Although subsequent tidal evolution 
tends to align the spin and the orbital angular momentum, the retrograde signature may still be 
perserved for relatively young binaries. 



6. Orbit and Spin Evolution of PSR J0045-7319/B-star Binary 



The PSR J0045-7319 binary, containing a 0.93 s radio pulsar and a massive B-star companion 
(M = 8.8M0, R = 6ARq) in an eccentric (e = 0.808) 51.17 days orbit (Kaspi et al. 1994), is 
unique and important in that it is one of the two binary pulsars discovered so far that have massive 
main-sequence star companions (The other one is PSR B1259-63 with a Be-star companion; 
Johnston et al. 1994). These systems evolve from MS-MS binaries when one of the stars explode in 
a supernova to form a neutron star. Thus the characteristics of such pulsar binaries can potentially 
be used to infer the physical conditions of neutron star at its formation. The PSR J0045-7319 
system, in particular, owing to its relatively small orbit and "clean" environment, exhibits two 
interesting dynamical orbital behaviors: (i) spin-orbit precessions due to the rapid, misaligned 
rotation of the B-star, which strongly suggests that the neutron star received a kick at birth from 
asymmetric supernova (Lai et al. 1995; Kaspi et al. 1996); (ii) Rapid orbital decay, on a timescale 
of -PorbZ-Porb = —0.5 Myr (shorter than the lifetime of the 8.8M0 B-star and the characteristic age 
of the pulsar) (Kaspi et al. 1996), which, together with a generic model of tidal evolution, can be 
used to constrain the age of the binary since the supernova and the initial spin of the pulsar (Lai 
1996). 

The mechanism for the orbital decay was discussed in Lai (1996). Since mass loss from the 
B-star is negligible (as inferred from dispersion measure variation; Kaspi et al. 1996b), the orbital 
decay must have a dynamical origin. The static tide, corresponding to the global, quadrupolar 
distortion of the star, cannot explain such rapid orbital decay because (i) it would require too 
short a viscous time (less than 30 days based on our numerical integrations of the tidal equations 
given by Alexander 1973) — Such a small viscous time is almost certainly impossible (e.g., even if 
the star were completely convective, the viscous time would still be longer than one year); (ii) It 
would lead to rapid spin-orbit synchronization and alignment even if viscosity were large enough 
to explain the observed -Porb) in contradiction with observations. It was suggested that dynamical 
tidal interaction between the pulsar and the B-star can explain the rapid orbital decay provided 
that the B-star has significant retrograde rotation with respect to the orbital motion so that the 
energy transfer is sufficiently large. The reason that retrograde rotation significantly increases 
the tidal strength is the same as discussed in §5.1, i.e., retrograde rotation shifts the resonance 
from high-order g-modes to lower-order ones, which couple more strongly to the tidal potential. 
Several simplifying assumptions were adopted in the preliminary analysis of Lai (1996). The 
most important one is that the mode damping times are assumed to be constant since only small 
numbers of modes are strongly excited. This assumption is likely to be inadequate, as correctly 
pointed out by KQ, who suggested that an additional ingredient, differential rotation, is needed to 
compensate for the longer damping times of lower-order g-modes. Our focus here, however, is still 
to explore the dyanmcal aspects of the tides more completely. 
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6.1. Dynamical Tidal Energy: Steady State 

The dynamical equations presented in §2 for the coupled orbital and mode evolution can be 
integrated numerically. Figure 7 depicts a typical example of such numerical integration. We start 
the calculation at the apastron (to = — Porb/2) with no initial oscillation, and include g-modes of 
order 2 — 8 (each has m = 0, ±2 components). The total angular momentum Jtrans and energy 
-E'trans transferred to the modes are calculated via 



trans 



* — GM'WimQa 
Jjl+l 

* — GM'WimQa . im<S> 



Jtrans(t) = / ^ " (-im) a^e^^ , (47) 



Obviously Jtrans is equal to the mode angular momentum J^, while -Etrans is equal to the mode 
energy Eg = X)q(1/2) [looP + ((^o — 2i?aCa)|aap] only in the absence of dissipation (r^ = 0). We 
have checked that angular momentum conservation ( + Jorb = const) and "energy conservation" 
(-2'trans + ^orb = const; of course the energy of the system, Eg + ii^orbi is not conserved due to 
dissipation) are satisfied numerically to high accuracy. 

Figure 7 reveals that, although initially the energy transfer during periastron passage varies 
from one orbit to another, after a few dissipation time, the mode energy {Ea) averaged over an 
orbit reaches a constant. In this steady state, the energy transferred to a mode at periastron is 
exactly balanced by the energy dissipated in one orbital period. Moreover, the mean mode energy 
{Ea) is approximately equal to the energy transfer in the first periastron passage, A£^q,, which has 
been discussed in §3. This apparent coincidence can be understood as follows. Consider a mode 
with a given a and sign of pattern speed. Its amplitude evolves according to the first term of 
Eq. (p7|), with the damping factor included: 



2z(a«-i3„)' Ju D^^\t') ■ 

Let tj = {2j — l)Porb/2 (with j = 0, 1, 2, • • •) be the times at apastron. After the fc-th periastron 
passage, the mode amplitude can be written as 

a„(tfe) = (Aac,)e(^'"--^")*^ y e(-^'^"+^")^'^-'^ = — — j-. i _ e(i<^..-r.)fePo.b (Aa„), (50) 

j=0 

where Aoq is the change of mode amplitude in the first periastron passage when there is no initial 
oscillation [cf. Eq. (p^)] 

GM'WimQa /-Wa ttGM'Q^ 

^'^^ = 1TU ^/ TmTi = -r^l+l, -Z-Klm[(7a)- (51) 

Clearly, when kPoj-y^Ta » 1, the amplitude aa{tk) becomes independent of k. The steady-state 
mode energy is given by 

^ ~ 2 [cOsh(rgPorb) - COs(a„Porb)] ~ 4 siu^ (cTaPorb/2) + (r„Porb)2 ' ^ ^ 



- 18 - 



where AEa = {GM''^R^ /Dp)TT^\QaKim{(Ta)\'^cra/{aa — Ba) is the energy transferred to the mode 
in the first periastron. Equation ( p2| ) estabhshes the relationship between steady-state tidal energy 
(Eg) = J2a{^a) and the initial energy transfer AEg = J2a^^a- Typically we have (Eg) ~ AEg, 
except when a mode is in close resonance with the orbit, i.e., = 27rA^/Porb = -^^^orb (where A'^ 
is an integer). 

6.2. Orbital Resonances 

We now consider the probability that the orbital evolution is driven by an near orbital 
resonance cXa-Porb = 2(A^ + £)tt, with |e| << 1. For the modes which contribute significantly to the 
energy transfer, cr„ ^ 217p, thus N = aa/^ovh ^ 2(1 + e)i/V(l - e)^''^ - 32. For (Ea) > WAE^ 
we require |e| < 0.05, and we use this to define the full width (APorbje (in the Porb-space) of 
the resonance: (APorb/-Porb)e = 2e/N. Suppose the energy transfer is dominated by a single 
mode a, then the time the orbit spends near resonance is {At)^ = (87r^e^/3A^)tdecay , where 
idecay = | -Porb /-forb | with Porb evaluated assuming {Eg) ~ {Ea) = AEa- On the other hand, the 
times it takes the orbit to evolve between resonances is (At),. = (2/A^)tdccay Clearly, the ratio 
{At)e/{At)r = (47r'^/3)e^ is much less than unity. Therefore it unlikely that the orbital decay is 
driven by such a close resonance (see also KQ). In the following, we shall set {Ea) — AE^- 

6.3. Retrograde vs. Prograde Rotations 

Figure 8 shows the dimensionless quantities [cf. Eqs. ( p^ and (p5|)] Tg^^ = T22^ (energy 
transfer in the rotating frame), 5*2 = S22 (angular momentum transfer) and T2 = T22 (energy 
transfer in the inertial frame) as a function of the stellar rotation rate (for f] = 7, corresponding to 
neutron star mass M' = 1.4 and orbital semimarjor axis 20R). All g-modes with n = 1, 2, • • • , 20 
are included in our calculations. As in prabolic tidal encounters (§5), we find that retrograde 
rotation increases the tidal energy significantly because energy transfer is shifted from higher-order 
to lower-order g-modes. We can see this effect clearly in Fig. 9, which shows contributions of 
different modes to the energy transfer for several different rotation rates: At Qg = 0, the dominant 
modes are ge-gs, while for Qg = —0.4, they are shifted to g3-g5. As a result, AEg for Ug = —0.4 is 
24 times larger than the nonrotating value. 

We see from Fig. 8 that moderate prograde rotation {0,g > 0) decreases AJg and AEg, and for 
sufficiently large Qg, both AJg and AEg can become negative, i.e., angular momentum and energy 
are transferred from the star to the orbit, leading to orbital expansion. We find that AJg < 
for ng > 0.319 and A^^ < for Qg > 0.362. These critical ilg's are larger than = 0.192, the 
orbital frequency at periastron 0. Although this result agrees qualitatively with that of KQ, there 



^The measured surface velocity of the B-star in the PSR J0045-7319 binary is 113 km/s (BeU et al. 1995). This 
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are significant quantitative differences: KQ found that AEg changes sign at 0^/0^ = 1.2 (from 
their Fig. 1 and and note that they used stehar radius R = QRq) whereas we find flg/^p = 1-9 0. 
Moreover, KQ found that at extreme rotation rate (say > 0.4) the magnitude of the energy 
transfer is larger than the nonrotating value (and they have concluded that prograde rotation can 
also explain the observed orbital evolution timescale, although with a wrong sign), whereas we find 
that |Ai?s| is much less than the nonrotating value even at ^Ig = 0.5. We think these differences 
mainly come from KQ's problematic treatment of mode decomposition (particularly the effect of 
Coriolis force). This is explained in the following paragraphs. 

KQ's results are based on the solution of an approximate mode amplitude equation in the 
rotating frame: 

(only I = 2 tides are considered). While they included the corrections to uja and Qa from the 
Coriolis force, KQ neglected the corresponding Coriolis force term, proportional to (d^ — iw^aQ,), in 
Eq. (53). Decomposing the external forcing 6"™*/!)'+^ into Fourier components with frequencies 



equal to multiples of Jlorbi we find that the steady-state mode amplitude is given by 



-imftst 



where Kim,k = ^orhKim{k^orh) ■ The mode amplitude in the inertial frame is simply e~ 
Note that Ki^^k is very small when m and k have the same signs, thus for m > we can neglect 
the positive-A: terms in the sum, and for m < neglect the negative-A: terms. Equation ( |54D is 
then equivalent to Eq. (1) of KQ (see also Kumar, Ao & Quataert 1995). However, it is easy to 
see that this expression is problematic: According to Eq. (^), even at the apocenter, the mode 
does not oscillate at its intrinsic frequency (cJq, in the rotating frame), and yet its amplitude still 
depends on the intrinsic mode quantities and Qa (which are dependent on the pattern speed 
of the mode). 

To see the problem more clearly, we can calculate the energy transfer and angular momentum 
transfer in the "first" periastron passage. As shown in §6.1, in the absence of close orbital 
resonance (§6.2), these are approximately equal to the steady-state mode energy and angular 
momentum. Applying the same procedure of §3 to Eq. (|53), we obtain the dimensionless tidal 



corresponds to a rotational angular frequency [in units of (GAf/i?^)^''^] of 0.22/ sinis„ (where isn is the angle between 
the spin axis and the line-of-sight). The component of fls perpendicular to the orbital plane is (0.22/ sinis„) cos/3, 
where the spin-orbit inclination angle /3 is constrained to be in the range 25° < jS < 55° or 125° < [3 < 155° (Kapsi 
et al. 1996; note that the other constraint, /3 < 41° or /3 > 139°, is effective only if we restrict the orbital precession 
phase angle to the first or the third quadrants). 

^This is different from static tide where Eorh and Jorh change sign at Sis — flp'. The weak friction theory (Alexander 
1973; Hut 1981) gives ^orb « -[1 - f{e)ns/np] and oc -[1 - g{e)ns/np], where f = 1, 0.979, 0.957 and 

g = 1, 1.216, 1.212 for e = 0, 0.808, 1. 
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energy and angular momentum: 



IT 



a 



a 



m 



(55) 
(56) 



T2 = Vr^y^— {0Ja-mQ.s)\K2rn{i^a-mQ.s)\'^ + {^a + mVLs)\K2-rn{i^a + rn9.s)? ■ (57) 



It is easy to see that Eqs. (|56D-(p7[) are identical to Eqs. (26) and ( p3| ) with Ba = —rnQg (i-e., 
with the Coriolis term neglected). For comparison, we also plotted in Fig. 8 the numerical results 
calculated using these expressions, which are in reasonable agreement with KQ's (see their Fig. 1). 
However, these expressions are certainly incorrect, as they would imply that a prograde mode and 
retrograde mode have the same frequencies in the rotating frame. Indeed, integrating Eq. (|5^ ) 
directly we find that after the first peristron passage, the mode amplitude (in the rotating frame) 
is a superposition of e«"^<^+*'^a* and g*™-''^"*'^"* — This is clearly unphysical. It is easy to show that 
the errorous expressions (55)-(|57|) lead to overestimates of the negative angular momentum and 



energy transfers for prograde stellar rotation. 



6.4. Long-Term Orbital Evolution 

In the steady state, the energy transferred to a stellar mode in each periastron passage is 
equal to the mode energy dissipated in one orbit, 2TaPoj.^{Ea) ~ 2raPovh^Ea. The corresponding 
angular momentum transfer is 2rQ,Porb^'^a- Thus the orbital energy and angular momentum 
decay rates are given by 



d, 

d{Jorh) V^OT- A T '^('^' 



^2r,(S„) ~ -^2r,A^„, (58) 
- -E2r„A.« = -^^ (59) 



If we define the mean damping time tdamp via J2a TaAE'ct = AEs/tdamp, then the current orbital 
decay rate is given by 

-(0.5Myr)-f:^)f^V (60) 



Porh ' ' VlO-V ltd 



amp 



Note that T2 = 10^^ corresponds to nearly maximum retrograde rotation. Thus to explain the 
observed -Porb; the damping times of the dominant modes (gs-gs; cf. Fig. 9) must be less than 
30 years. As pointed out by KQ, such a short radiative damping time is unlikely if the star is 
rigidly rotating, but would be possible if there is differential rotation. Such differential rotation is 
naturally expected since tidal torque deposits angular momentum mainly in the region near the 
stellar surface (Goldreich &; Nicholson 1989). 
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Independent of the mode damping rates, equations ( |5q ) and (|59D provide a scaling relation for 
the orbital decay rate as a function of the orbital and spin parameters. For 4 < r/ < 10, we can fit 
the function T2 by T2 oc r]~'^'^, where v depends mainly on Qg but only slightly on the eccentricity 
(as long as e ^ 0.7): we find = 1, 0.5, 0.2 for fig = 0, -0.2, -0.4 (see Fig. 4). We then have 



-forb 
-forb 
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(61) 



where the proportional constant can be fixed by the current observed value of Porb- Similar 
equations for e and fig are easily obtained. These equations can be used to infer the long-term 
evolution of the binary system. Note that since A^;^ ~ P{GM/R^y/^AJg, where /? is of order 
unity (for the current system, 77 = 7, we find /3 = 3 for = and /? = 5 for Qg = —0.5), the rates 
at which Jorb and Jg change are given by 



3 



1 



1 + e 



nr. 



-forb 




Js 


0.50/3(1 - e) 


-forb 


5 


Js 


A(l + e)V3 



orb 



-forb 



(62) 



where we have taken the moment of inertia / = O.IXMR^ , and and Jlp are in units of 
{GM/R^y/^. Both I Jorb/^rbl and \jg/Jg\ are much smaller than l-PorbZ-Porbli i-e., the timescale for 
changing the orbital angular momentum and the timescale for synchonizing and aligning the stellar 
spin are much longer than the orbital decay time Q The orbit therefore evolves with decreasing 
orbital period and eccentricity, while the periastron distance Dp remains approximately constant. 
Since AEg mainly depends on Dp, the orbital decay rate l-Porb/^orbl oc AE's/|£'orb| was larger at 
earlier times (when l-Borbl was smaller). This implies that there is an upper limit to the age of 
the binary since the neutron star formation. Integrating the parametrized evolutionary equations 
backward in time, we find that regardless of the various uncertainties, the orbital age is always 
less than 1.4 Myr (Lai 1996). This is significantly smaller than the characteristic age (3 Myr) of 
the pulsar, implying that the latter is not a good age indicator. The most likely explanation for 
this discrepancy is that the initial spin period of the pulsar is close to its current value. Thus the 
pulsar was either formed rotating very slowly, or has suffered spin-down due to accretion in the 
first ~ 10^ years (the Kelvin-Holmholtz time of the B-star) after the supernova (E. van den Heuvel 
1996, private communication). 



7. Discussion 

The main result of this paper has been summarized in the abstract. Here we simply 
note that our study of dynamical tides in rotating stars is based on an approximate scheme 



^This is different from evolution driven by static tide, for which -Eorb — ilpJorb, and we can show that in this 
case I jorb/Jorb| and |Ja/Js| are larger than those given in Eq. ( |62| ) by a factor f2p ^ (although both are still less than 

|Porb/Porb|). 
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of mode decomposition, which we beheve gives physically meaningful results. Nevertheless, a 
mathematically rigorous theory describing the response of individul modes in rotating stars to 
external forcing is highly desirable, and the role of the continuous mode spectrum (Dyson & 
Schutz 1979) need to be clarified. 

In addition to the two applications (tidal capture binary formation and the orbital evolution of 
pulsar binaries) examined in this paper, tidal excitation of oscillation modes may also be important 
in the final phase of inspiraling neutron star binaries. Extracting gravitational wave signals from 
noise requires accurate theoretical waveforms in the frequency range 10 — 1000 Hz, corresponding 
to the last few minutes of the binaries' life (Cutler et al. 1993). For nonrotating neutron star, 
the orbital phase error induced by resonant excitaion of g-mode is negligible (Reisenegger & 
Goldreich 1994; Lai 1994). With rapid rotation, the f-mode frequency in the inertial frame can 
be significantly reduced, therefore resonant excitation of f-modc is possible even at large orbital 
separation. We wish to study this and related issues in a future paper. 

I thank Pawan Kumar for sending me his papers prior to publication and several interesting 
conversations. I also thank Lars Bildsten, Lee Linblom and Greg Ushomirsky for discussing modes 
in rotating stars. This research is supported by the Richard C. Tolman Fellowship at Caltech, 
NASA Grant NAG 5-2756, and NSF Grant AST-9417371. 



A. Equations for General Spin-Orbit Inclination Angle 

In the main body of our paper we assume for simplicity that fig is aligned with the orbital 
angular momentum. Generalization to the cases of arbitrary spin-orbit inclination angle is 
straightforward. Define a coordinate system (xyz) centered on M with the z-axis along the orbital 
angular momentum and the x-axis pointing to the pericenter in the orbital plane. Define another 
coordinate system {x'y'z') with the z'-axis along ftg and the y'-axis in the orbital plane. Let the 
angle between the z-axis and ^'-axis be /3 (the spin-orbit inclination angle), and the angle between 
the y-axis and y'-axis be a. Thus the {x'y'z') frame is related to the (xyz) frame by Euler angle 
(a,/3,7 = 0). The function Yijn{6,(p) is related to Yi^> {6' , (p') by 

YlmiOA) = J^^^Xi^^mm'ie', <!>'), (Al) 
m' 

where the Wigner P-function is given by 

P«^(«,/3) = e^"^'"[(Z + m)!(/-m)!(Z + m')!(/-m')!]'/' 



k 



^ W _l_ m _ bW (I -L^' — (h 
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(e.g., Wybourne 1974). The mode amplitude equation (^) should then be replaced by 

■• I 2 o-D /• ■ ^ I o • GM'Wim'Qal^{l) -im'<S> / a q\ 

tta + (T^aa - 2tBa{aa - tCTaaa) + ^laaa = ^-mm'^ ■ (^^j 



It is easy to show that in the egneral cases, equations (^)-(|3^) still apply, provided that Kimica) 
is replaced by 

m' 

[recall that a = {njm}]. Note that in general, the m = ±1 modes also contribute to the tidal 
excitation (unlike the aligned case, where only the m = 0, ±2 modes are excited). 



B. The Functions Qai and Ba 

With ^r(r) and 5P{r) expressed in Eq. (^), the other eigenfunctions are 

6p{r) = 6p{r)He''^'t' = ( ^ + ^^r] He'"'^ (Bl) 

\ 9 J 

Ur) = -^('4 + ^.H)e^-t (B2) 



1 — q'^fi'^ V d6 sin( 

1 — V o6 smO J 

where the mode index a = {njm} has been surpressed, q = 2Qs/L0a, ^ = 0086*, and 
^^(r) = {6P)/{pruj^). We normalize the Hough function via JdnH^{9) = 2tt JdfiH'^ = 1. The 
eigenfunction is normalized according to 

I d^xpi* ■ i= dr r^p + A^i) = 1, (B4) 



(B5) 



where 

Note that A— >/(/ + !) as >0. The tidal coupling coefficient defined by Eq. (^) is 

Qci = Qe j^drr'+Hp^{r), Qe = j di}He-'""^Yi^. (B6) 
The function defined by Eq. (^) is 

Ba = — mQs + 2i^ls J d^x pS,'^{(,e COS 6 + sin 6) 

= -mVLs + 2ns (^Bei j\r pi\ + Bq^ £ dr r^p^rU^ , (B7) 
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where 



Be2 



2tt 
2tt 



-i(l-(7V)' 



qfi{l - I?) 



H - m{l + q jJ. )H—— 



\dn J I- fi 



+ 



.2^2) 



9iJ 



mi7^ - qa(l - U^)H— 
on 



Note that Bgi m/2 and Bgi ^ m as q ^ 0. 



(B8) 
(B9) 



C. Asymptotic Expressions for Kim{io) 

For parabohc orbit and in the hmit y = -q a » 1 (where a = u) — rnQg is the mode frequency 
in the inertial frame), the functions Kim{cr) can be evaluated analytically. Following PT, we write 



V2Wim . , . 



IT 



Iim{y)= / d2;(l + x^)"'cos [2^/^y(x + x^/3) + 2mtan~^; 
Jo L 



(CI) 



(C2) 



Using the method of steepest decents, the integral ( p2D can be evaluated in the asymptotic limit 
{y » 1). For 1 = 2, the leading terms are 



ho{y) 



^z3/2exp(-2z/3) { 1 - ^z-V2 + 



-^z I exp( 
7rV2 



4 



-2^/3) (^1+ 2 



32 



-1/2 



exp(-2z/3) 1 



48 



(C3) 
(C4) 
(C5) 



where z = \f2y and we have assumed y > 0. Comparison with numerical integration indicates 
that these asymptotic expressions are accurate to within 2% for y > 2. 
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Figure 1 - Frequencies uj (in the rotating frame) of selected j = 2 (correcponding to / = 2 in 
the — > limit) g- modes as a function of rotation rate $7^. Both tu and 0,s are expressed in 
units of u!^^\ the corresponding mode frequency at zero rotation. The solid curves are for m = 0, 
short-dashed curves for m = —2, and long-dashed curves for m = 2. For each m, the four curves 
correspond g-modes of different radial order: g2, gs, gio [these curves terminate at J7s = 0.5], and 
goo (the WKB limit). The numerical results shown in this figure can be fitted by Eqs. (^)-(|4l]). 

Figure 2 - The tidal coupling coefficients Qa = Qa2 of selected j = 2 modes as a function of 
rotation rate ^Ig (in units of uj^^\ the corresponding mode frequency at zero rotation). The solid 
curves are for m = 0, short-dashed curves for m = —2, and long-dashed curves for m = 2. For each 
m, the three curves correspond g-modes of different radial order: gs, gio (these curves terminate 
at = 0-5) and g2o (Note that the gs curves almost overlap with the gio curves). 

Figure 3 - The frequencies lo and angular coupling coefficients Qe of selected j = 4: g-modes 
(heavy curves) compared with those of the j = 2 modes (lighter curves). The solid lines are for 
m = 0, short-dashed lines for m = —2 and long-dashed lines for m = 2. (a) gives the frequency (in 
the rotating frame) oj in the WKB limit {n — > oo), and (b) shows the quantity Qe as defined in 
Appendix B. 

(r) (r) 

Figure 4 - Dimensionless functions T2 , ^2 and T2 = T2 + ^382 [where fls is the rotation 
rate in units of (GM/R^)^^"^] as a function of i] for = (solid lines), 0.2 (dotted lines), —0.2 
(short-dashed lines), 0.4 (long-dashed lines) and —0.4 (dot-dashed lines) in a parabolic encounter. 
The star is modeled as a P = 4/3, Fi = 5/3 polytrope. Note that for = 0.4 (long-dashed lines), 
S2 becomes negative for t] > 5.36 [see (b)], and T2 becomes negative for r/ > 6.45 [see (c)]. 

Figure 5 - Energy transfer as a function of periastron distance Dp during an encounter between 
a main-sequence star M = IMQ, R = IRq (modelled as a F = 4/3 polytrope) and a M' = 1.4M0 
point mass (neutron star). The four horizontal lines show the kinetic energy at infinity fj,v'^/2 of 
the two stars with relative velocity v = 20, 10, 5, 2.5 kms^^. The solid curve is for Qg = 0, the 
dotted curve = 0.2, the long-dashed curve Qg = 0.5 (close to the maximum rotation rate), the 
short-dashed curve ^Ig = —0.2, and the dot-dashed curve ^Ig = —0.5. 

Figure 6 - Energy transfer as a function of periastron distance Dp during an encounter between 
a main-sequence star M = IMQ, R = IRq (modeled as a F = 5/3 polytrope) and a M' = 1.4Mq 
point mass (neutron star). The four horizontal lines show the kinetic energy at infinity fj,v'^/2 of 
the two stars with relative velocity v = 20, 10, 5, 2.5 kms~^. The solid curve is for Qg = 0, the 
dotted curve 0,g = 0.3, the long-dashed curve 0,g = 0.6 (close to the maximum rotation rate), the 
short-dashed curve ^Ig = —0.3, and the dot-dashed curve ^Ig = —0.6. 
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Figure 7 - Evolution of (a) the energies SEorh = -E'orb — Eo (wliere Eq is the initial orbital energy) 
and E^"^^ (the mode energy in the rotating frame), and (b) the angular momenta 5Jorb = >^rb — Jo 
(where Jq is the initial orbital angular momentum) and SJg (the mode angular momentum) 
due to dynamical tide with = —0.4. Only g2-g8 modes are included in the calculation, and 
relatively large damping rates are chosen for clearer illustration: Tq, = 0.1/Porb (heavy lines) and 
r„ = 0.05/Porb (light lines). 

Figure 8 - Dimensionless tidal energies Tg*^^ (in the rotating frame), T2 (in the inertial frame) 
and angular momentum 5*2 as a function of the stellar rotation rate ^Ig [in units of (GAf/i?'^)^/^]. 
The dimensionless orbital parameters are r] = 7 and e = 0.808. The solid curves are obtained 
using the correct expressions (P)-(|3|); the dashed curves are obtained using the same equations 
but assuming = —mil.g (i.e., neglecting the Coriolis term) — These demonstrate that our 
results are insensitive to uncertainties in Eqs. (|3ll)-(^). Note that S2 and T2 become negative 
for 0,s > 0.311 and ^Ig > 0.356 respectively. The dotted curves are obtained using the incorrect 
expressions (|55D-(|57|); these give results similar to those of KQ. 

Figure 9 - The contributions of g-modes of different order (n) to the function T2 [the dimensionless 
energy transfer in the inertial frame; cf. Eqs. (p2D, (|3l|)] for 17^ = (triangles), = —0.4 (squares) 
and Qs = 0.4. The solids circles indicate positive contributions, the open circles indicate negative 
ones. 
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